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ABSTRACT 

We analyze a sample of galaxies with stellar masses greater than IO'^Mq and with redshifts in 
the range 0.025 < z < 0.05 for which HI mass measurements are available from the GALEX 
Arecibo SDSS Survey (GASS) or from the Arecibo Legacy Fast ALFA survey (ALFALFA). 
At a given value of M», our sample consists primarily of galaxies that are more Hl-rich than 
average. We constructed a series of three control samples for comparison with these Hl-rich 
galaxies: one sample is matched in stellar mass and redshift (C^*), the second sample is 
matched in stellar mass, NUV - r colour and redshift {CM*.NUv-r), and the third sample is 
matched in stellar mass, NUV - r colour, stellar surface mass density /j, and redshift. We 
generated self-consistent 7-band photometry (FUV, NUV, u, g, r, /, z) for all galaxies and we 
used this to derive inner colours, outer colours, asymmetry and smoothness parameters. We 
also used standard SED fitting techniques to derive inner and outer specific star formation 
rates. As expected, Hl-rich galaxies differ strongly from galaxies of same stellar mass that 
are selected without regard to HI content. The majority of these differences are attributable 
to the fact that galaxies with more gas are bluer and more actively star-forming. In order 
to identify those galaxy properties that are causally connected with HI content, we compare 
results derived for the HI sample with those derived for galaxies matched in stellar mass, size 
and NUV-r colour The only photometric property that is clearly attributable to increasing HI 
content, is the colour gradient of the galaxy. Galaxies with larger HI fractions have bluer, more 
actively star-forming outer disks compared to the inner part of the galaxy. Hl-rich galaxies also 
have larger §-band radii compared to /-band radii. Our results are consistent with the "inside- 
out" picture of disk galaxy formation, which has commonly served as a basis for semi-analytic 
models of the formation of disks in the context of Cold Dark Matter cosmologies. The lack 
of any intrinsic connection between HI fraction and galaxy asymmetry suggests that gas is 
accreted smoothly onto the outer disk. 
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1 INTRODUCTION 

The fraction of the baryons in dark matter halos that are locked up 
in stars in the central galaxy reaches a maximum of 20%, at masses 
somewhat below that of the Milky Way, and falls rapidly at both 
higher and lower masses (Guo & White 2010). In the absence of 
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feedback effects from supernova or AGN , most of these baryons are 
expected to cool, accrete and form stars jKaufFmann et al.|[l993l) . 

There is some observational evidence that gas from the exter- 
nal environment continues to accrete onto galaxies at the present 
day. Without a continuous supply of gas, most star fo rming galax- 
ies w ould run out of gas on timescales of a few Gyr jLarson et al.l 
Il980l) . HI cloud complexes, Hl-rich dwarfs in the vicinity of spiral 
galaxies, extended and warped outer layers of HI in spiral galaxies, 
and lopsided galaxy disks have all been cited as evidence for ongo- 
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ing gas accretion in nearby galaxies jSancisi et al .I2OO8I) . However, 
estimates of the accretion rates in the form of neutral hydrogen in 
nearby spirals consistently give values much lower than those re- 
quired to sustain star formation at their observed rates (see Frater- 
nali (2010) for a recent review). This leads to the hypothesis that 
most of the baryons that reside outside galaxies at the present are 
in the form of ionized gas (e.g. Fukugita, Hogan & Peebles 1998). 
There is also some evidence from observations of ionized silicon in 
high and intermediate-velocity clouds in the Milky Way, that gas in 
a low-metallicity ionized phase in the halo can provide a substan- 
tial ( 1 Mo/yr) cooling inflow to replenish star formation in the disk 
(Shull et al. 2009). 

One interesting question that has not been considered very 
much in the literature, is whether gas accretes onto all galaxies in a 
smooth, continuous fashion, or whether accretion is more stochas- 
tic. If accretion has been stochastic, then galaxies that have recently 
acquired gas from the external environment might be expected to 
have unusually high HI mass fractions. The neutral gas content of 
a galaxy has long been known to correlate with its other physi- 
cal properties. HI mass fractions increase smoothly along the Hub- 
ble s equence from the early-type (SO) to the late-type (Im) end 
( [Roberts & Haynes 1994). There are also correlations between HI 
mass fraction and galaxy properties such as stellar mass (M,), stel- 
lar mass surface density (/j.) and broadband colour. These corre- 
lations form the basis for the so-called "photometric g as fraction" 
technique for predicting the HI gas content galaxies jKannappaiil 
[2004 ; Zhang et al. 2009). Rather little attention has been given to 
the scatter in gas fraction at a given value of M,, ji, or colour. Re- 
cently, Zhang et al (2009) found that more gas-rich galaxies were 
systematically deficient in metals at fixed stellar mass and they in- 
terpreted this as evidence for recent cosmological infall of gas in 
these systems. However, up to now there has been no systematic 
study of how galaxy properties vary as a function of HI over- or 
under-abundance. 

One problem that has hampered progress on this front has 
been the lack of suitable data sets. Blind HI surveys, such as the HI 
Parkes All Sky Survey (HIPASS, Barnes et al 2001) or the Arecibo 
Legacy Fast ALFA survey (ALFALFA, Giovanelli et al 2005) have 
produced large, unbiased samples of galaxies selected by HI mass. 
However, these surveys are still relatively shallow, so the fraction 
of Hl-poor galaxies is small. 

Recently, Catinella et al (2010, hereafter CIO) explored how 
HI gas fraction scales as a function of galaxy stellar mass, galaxy 
structural parameters and NUV-r colour, using data from the 
GALEX Arecibo SDSS Survey (GASS). The survey is an ongo- 
ing large programme at the Arecibo radio telescope that is gather- 
ing high quality Hl-line spectra for an unbiased sample of ~ 1000 
galaxies with stellar masses greater than IO'^Mq and redshifts in 
the range 0.02 5< z < 0.05, sele cted from the Sloan Digital Sky 
Survey (SDSS. lYork eTal] (I2OO0I)) spect r oscop ic and Galaxy Evo- 
lution Explorer (GALEX, iMartin et al.l ( l2005h ) imaging surveys. 
The galaxies are observed until detected or until a low gas mass 
fraction limit (1.5-5 % ) is reached. CIO quantify in detail how the 
mean atomic gas fraction of the galaxies in their sample decreases 
as a function of stellar mass, stellar mass surface density, galaxy 
bulge-to-disc ratio (as measured by the concentration index of the r- 
band light), and global NUV-r colour. The fraction of galaxies with 
significant (more than a few percent) HI decreases sharply above a 
characteristic stellar surface mass density of 10**^ M© kpc"-. The 
fraction of gas-rich galaxies decreases much more smoothly with 
stellar mass. 

In this paper, we extend the study of CIO to study how the 



resolved UV/optical photometric properties of galaxies depend on 
atomic gas fraction. We have developed a photometric pipeline 
that transforms the SDSS and GALEX images to the same geome- 
try and effective resolution, so that our photometric measurements 
trace the same part of each galaxy at different wavelengths (Wang 
et al 2009). We use these images to study UV-optical colours in 
their inner and outer regions of the galaxies in our sample and to 
study how colour (and hence specific star formation rate) gradients 
depend on HI content. In order to extend the range in M(HI)/M, 
that we are able to probe, we combine the sample analyzed in CIO 
with galaxies from the GASS "parent sample" (described in Sec- 
tion O which have HI detections in the ALFALFA survey. This 
considerably improves the statistics for galaxies with large values 
of M(HI)/M,, which are relatively rare. In order to isolate the effect 
of the presence or absence of gas on the photometric properties of 
galaxies, we compare our results with three control samples, which 
have similar optical properties to the galaxies in our sample with HI 
mass measurements, but are selected without regard to HI content. 
These control samples are described in detail in Section 1231 Sec- 
tions 13. H and I3.3| present the methods used to measure our photo- 
metric parameters, and in Section 4 we discuss the trends in galaxy 
size, NUV-r and specific star formation rate gradients, and higher 
order morphological parameters such as asymmetry and smooth- 
ness as a function of HI mass fraction. Finally Section [5] summa- 
rizes our results. Throughout this paper, we assum e a cosmology 
with //o=70 km s"' Mpc"', U„=0.3, and aA=0.7 jTegmark et aU 
I2OO4) . A Chabrier initial mass function is assumed in the stellar 
population synthesis analysis ( IChabrierll2003l) . 



2 DATA 

2.1 The parent sample 

We select galaxies with stellar masses M, > 10'°Mo in the redshift 
range 0.025 < z < 0.05 from the sixth data release (DR6) of the 
SDSS Survey, which lie within the maximal ALFALFA footprint. 
We match this catalogue with the fourth data release (GR4) of the 
GALEX survey and this yields a sample of 10468 galaxies. We 
further limit the sample to face-on galaxies with h/a > 0.4 (or an 
inclination of less than 66.4 deg) where b and a are the minor and 
major axes of an ellipsoidal fit (from SDSS r band) to each galaxy. 
This yields a final sample of 8429 galaxies, which we refer to as 
the parent sample from now on. We restrict the sample to face- 
on systems, so that our estimates of star formation rate, which are 
derived from the UV/optical photometry, will be less affected by 
dust attenuation effects (Section[33). 



2.2 The sample of galaxies with HI mass measurements 

This sample consists of galaxies in the parent sample for which 
we have catalogued HI detections from the 40% ALFALFA sur- 
vey (q'.40, which w ill be discussed in iMartin et al.l ( I2OI0I) and 
iHavnes et aTl ( l2010h ). as well as galaxies that are included in the 
first data release (DRl) of the GASS survey (see CIO for a de- 
tailed description). Our HI sample consists of 458 galaxies with 
a median HI mass fraction M(HI)/M, of ~ 30%. Figure [T] shows 
the distribution of M(HI)/M, for these galaxies. The left-top panel 
of Figure |2] show a gallery of g, r, / colour images generated using 
the SDSS Image List ToolQ for galaxies with stellar masses in the 
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range 10'"'' to 10" ' Mq and HI mass fractions (log M(HI)/M.) in 
the range -0.84 to -0.48 dex. We note that the average value of log 
M(HI)/M, for galaxies in this stellar mass range is — 1.4 (see CIO), 
so these galaxies are all significantly more Hl-rich than is typical. 

2.3 Control samples 

We limit the galaxy pool from which we draw the control galaxies 
to the regions of sky not yet covered by our current ALFALFA 
or GASS catalogues. The sample from which we draw the control 
galaxies comprises 5527 galaxies from the "parent sample". We 
create three control galaxy samples with similar optical properties 
to the galaxies with HI detections. 

The first control sample {Cm,) is matched only in redshift and 
stellar mass. For each galaxy with an HI detection, we search within 
a distance of 0.028 dex in M, and 0.005 in z for the nearest neighbor 
in the plane of stellar mass versus redshift. Each galaxy only enters 
the control sample once. We repeat this procedure twice, increasing 
the AlogM, tolerance to 0.03 and 0.035, so that in the end, each 
galaxy with an Hl-detection has 3 galaxies matched in stellar mass 
and redshift. The bottom-left panel of Figure[2]shows a montage of 
Cm, galaxies that have been matched to galaxies in Rows 1-2. Note 
that because the control galaxies are matched in redshift, we can 
meaningfully compare the morphologies, sizes and colours of the 
HI rich and corresponding Cm, galaxies in this figure. 

Figure [3] compares the distributions of redshifts (z), stellar 
masses (log M,), total NUV - r colours, half-light radius (i?5o), 
local environment densities Q (p), concentration indices (Rgo/Rso), 
g — i colour difference (A(g - /), see Sect.3.1) and 25 mag arcsec"- 
isophote diameters (D25) for galaxies with HI measurements (open 
histograms) and for Cm, galaxies (hatched, red histograms). The 
logarithm of the KS-test probabilities that both histograms are 
drawn from the same underlying distribution are noted in the cor- 
ner of each panel: a value of -2 means there is 99% confidence that 
the null hypothesis that the two samples are drawn from the same 
distribution can be rejected. 

Figure[3]confirms the visual impression gained from compar- 
ing the left-top and left-bottom panels of Figures |2l Galaxies in 
our sample with measured HI masses are, on average, significantly 
more gas rich than Cm, galaxies and are also found to have bluer 
global colours, lower concentrations, larger sizes, and stronger neg- 
ative colour differences (i.e. they are bluer on the outside). Galaxies 
with HI measurements are found in lower density environments, 
which is not too surprising, since galaxy colour and environment 
are known to be strongly correlated. 

We have constructed a second control sample, CM,,Nuv-r-, that 
is matched in NUV - r colour in addition to stellar mass and red- 
shift. For each galaxy with an HI measurement, we search within 
a distance of 0.18 dex in log M», 0.31 mag in NUV - r colour 
and 0.005 in z for the nearest neighbor in the stellar mass-redshift- 
colour plane. Note that the chosen tolerances in log M, and NUV-r 
colour are comparable to the errors in these quantities. We repeat 
the above process again, with the log M, distance increased to 0.25 
and the (NUV - r) distance increased to 0.43 mag. In the end 
each HI galaxy has 2 control galaxies matched in stellar mass and 
(NUV - r) colour. 



Note that p is calculated as the number of tracer galaxies inside a cylinder 
of diameter 2 Mpc h"' and length 16 Mpc h"' divided by the mean number 
of galaxies inside such a cylinder at that redshift. The p values have been 
extracted from the NYU Value-Added Galaxy Catalog jfilanton et al.l2005l) 



The CM,,Nuv-r control sample (examples are shown in the 
right-top panel of Figure |2j allows us to investigate the extent to 
which the NUV - r colour can serve as a "proxy" for the HI con- 
tent of a galaxy. Note that this is the underlying assumption of many 
photometric gas fraction techniques - if HI fraction and colour 
are exactly equivalent, then the HI sample and the CM,,Nuv-r con- 
trol sample should have identical properties. Figure |4]presents his- 
tograms of the properties of these two samples. As can be seen, 
there is no longer any significant difference in concentration in- 
dex between the sample with HI masses and the CM,,Nuv-r control 
sample, but the size differences do persist. The differences in the 
distribution of local environment parameters and colour gradients 
decrease significantly, but are still significant. 

Finally, we have constructed a third control sample, 
C M,.Nuv-r.ii, > that is matched in half-light radius measured from the 
SDSS !-band image (RsoU)) in addition to stellar mass, NUV - r 
colour and redshift. For each galaxy with an HI measurement, we 
search within a distance of 0.3 dex in log M,, 0.52 mag in NUV-r 
colour, 1.9 kpc in Rsoii), and 0.005 in z for the nearest neighbor 
in the stellar mass-redshift-colour-size plane. In this case, there is 
only one control galaxy matched to each HI galaxy. Catinella et al 
(2010) showed that a HI mass could be most accurately "predicted" 
using a linear combination of NUV-r colour and log^,. Figure 5 
shows that very few apparent differences between the HI sample 
and the the CM,.Nuv-r,ii, control sample now remain. There is no 
longer any significant difference in the distributions of local envi- 
ronment parameters. There does appear to be a difference between 
the distribution of g - i colour gradients, which we will investigate 
in more detail in the next sectior0. 



3 MEASUREMENTS 
3.1 2-zone photometry 

In this section, we analyze how the sizes and colour gradients of 
galaxies of fixed stellar mass depend on atomic gas fraction. We 
have generated two different sets of aperture-matched photometry 
for galaxies with HI measurements and also for the galaxies in our 
control samples. The first set of measurements is generated from 
GALEX FUV and NUV images, and SDSS images in the u, g, r, i 
and z-bands. Images in the FUV, g, i, and z bands are regis- 
tered so that they have the same geometry as the NUV image, and 
are then convolved to the same effective resolution as the NUV 
image. All neighboring objects are masked and the r-band image 
is u sed as the referenc e for flux measurements in all other bands 
(see I Wang et"ai] | |2009|) for details). We use this set of photometric 
measurements in our analysis of UV/optical colours and colour gra- 
dients, and for self-consistent SED fitting using all 7 photometric 
bands (FUV,NUV, u, g, r, /, z). 

The second set of measurements is generated from the SDSS g 
and / band images, which have significantly higher resolution than 
the GALEX images (~ 1 " instead of 5 "). The /-band image is 
registered to match the geometry of the g band image, and the two 
images are convolved so that they match to a common PSF. The 
g-band image is used as the reference image for flux measurements 
in both bands. This set of photometric measurements is used when- 
ever we consider only optical colours, sizes or colour gradients. 

^ We also made a control sample similar to CM,.NUV-r, but which ex- 
cl uded galaxies hosting a ctive or composite nuclei (AGN) using the criteria 
in lKauffman et aljj2003l) . We did not find any difference in the results. 
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Figure 1. The distribution of HI mass fractions for the 458 galaxies in our sample with HI mass measurements. 
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Figure 6. The convolved r-band image of the GASS galaxy G3759. The 
dashed line shows the half-light ellipse, the dashed-dotted line shows the 25 
mag/arcsec^ isophote, and the solid line shows the ellipse at 2.5 times the 
Kron radius. The region within the dashed line is defined to be the "inner" 
region, while the region between the dashed and the solid lines is defined as 
the "outer" region. 



We divide each galaxy into an inner and an outer region: the 
former is enclosed by R^q, the radius enclosing half the total r-band 
light (determined from the PSF-convolved r band image), and the 
latter is defin e d as t he region between i?5o and 2.5 times the Kron 
radius jKrorj ( Il980h . also determined from the PSF-convolved r 
band image). We measure fluxes in all the bands for both regions. 
We define the (NUV- r) colour difference A(NUV- r) as (NUV - 
r)mit - (NUV - r)i„ (likewise for A(g - /)), so that negative values of 
A imply that the outer region of the galaxy is bluer than the inner 
regior|j. 

To illustrate. Figure [6] shows the Hl-detected GASS galaxy 
G3759. The dashed curve indicates the ellipse enclosing half the 
r-band light, the dot-dashed curve indicates the 25 mag/arcsec- 
isophote, and the solid curve the ellipse with a semi-major axis 
equal to 2.5 times the Kron radius. 



' We have checked that if we normalize the colour difference A(NUV — r) 
by the size of the galaxy, the results presented in this paper do not change. 



3.2 Image stacking 

In Appendix [A] we demonstrate that the 2-zone (NUV - r) color 
measurement is reliable only when the radius of the inner zone is 
larger than 6 arcsec (about the FWHM of the PSF of a GALEX 
NUV image) and the error on the total (NUV-r) colour (a(NUV - 
r)) is smaller than 0.1 mag. If we were to apply both a size and an 
error cut, we would be left with ~ 74% of the galaxies from the 
HI sample, ~ 36% of the galaxies from the Cm* sample, ~ 70% of 
the galaxies from the Cm..jv(/v-, sample and ~ 74% of the galaxies 
from the CM,,Nuv-r.nt sample. If we only apply the cut on size, 78% 
of the HI sample, 58% of the Cm. sample, 76% of the Cut-Nuv-r 
sample and 80% of the CM..jv(/v-r.^. sample remain. 

One problem with making a cut in a(NUV - r) is that the re- 
maining galaxies are biased toward NUV-bright (i.e. blue) objects, 
especially in the Cm* sample. In order to avoid this problem, we cut 
the sample only by size, and we stack the images in bins of stellar 
mass and HI gas fraction in order to boost the S IN of the colour 
measurements. We then measure average colour and SFR gradients 
for the stacked images. The increased 5 IN of the stacked images 
will also improve the quality of the SED fitting (Section[33](. 

We first select galaxies which have Rsa larger than 6", from the 
HI and "parent" samples. We then make new C^,, CM,,Nuv-r and 
CMt.Nuv-r,iit control samples as described in Section [231 Because 
the "galaxy pool" from which we draw these control galaxies is 
now smaller, each galaxy from the HI sample only has one con- 
trol galaxy matched in stellar mass and one control galaxy matched 
in stellar mass and (NUV - r) colour. For simplicity, we will still 
call these three new samples the HI sample. Cm, sample, CMt.Nuv-r 
sample and CM*.Nuv-r.ii* sample. 

We divide the new HI sample into 4 stellar mass bins with log 
M./Mo ranges of 10-10.25, 10.25-10.5, 10.5-10.75 and 10.75-11.5. 
Each stellar mass bin includes 80-90 galaxies. In each stellar mass 
bin, we further divide the galaxies into two groups at the median 
value of M(HI)/M,. In total, the HI sample is divided into 8 groups 
for stacking, with 40-50 galaxies in each group. We also divide the 
matched Cm. and CMt.Nuv-r samples into 8 groups for stacking; 
each group in these control samples will then be matched to the 
corresponding group from the HI sample. 

Our stacking procedure consists of the following steps. We 
rescale the images to the median value of R;,q for the galaxies in the 
group. We subtract the background, correct for Galactic extinction 
and any offset in photometric zeropoint, align the centers of all the 
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Figure 2. A montage of g, r, i colour images of galaxies selected from the HI sample (top-left panel), and the corresponding control galaxies from the Cm, 
(bottom-left panel), CMt,NUV-r (top-left panel) and CM-,NUV-r,ii. (bottom-right panel) samples. The scale of each image is ~0.2" pixel"'. The galaxies are 
selected to have stellar masses in the range from lO'""^^ to lO""' M©. The galaxies from the HI sample are ordered by increasing logM(HI)/M,, from -0.84 
to -0.48 
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Figure 3. Comparison of distributions of redshifts (7). stellar masses (log M,), total NUV — r colours, half-light radii (R50), local environment densities (p), 
concentration indices (Rgo/Rso)^ S ~ ' colour differences {A(g - i)) and 25 mag arcsec"^ isophote diameters (Dij) for the sample with HI measurements (open 
histograms) and the Cm* control sample (red, hatched histograms). The logarithm of KS-test probabilities that the two histograms are drawn from the same 
underlying distribution are denoted in the comer of each panel. 




Figure 4. Same as Figure|3] but for the CMt,NUV-r control sample. 




Figure 5. Same as Figure|3] but for the CMt.NUV-r,iit control sample. 
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images, and then add them together. We create stacked images for 
the images (convolved to the resolution of the NUV image) for all 
7 photometric bands. The stacked r-band image is used as the ref- 
erence image to measure photometric parameters. We note that the 
stacked images all have a (t(NUV - r) less than 0.01. 

3.3 2-zone specific star formation rates 

We have used a spectral energy distribution (SED) fitting technique 
to derive the average specific star formation rates (sSFR) in the in- 
ner and outer reg ions of our stacked sample of galaxies. We follow 
the method from lSalim et al. (2007), hereafter S07, with a few key 
chang es that are detailed below. S07 used the iBruzual & Charloll 
( l2003h population synthesis code to create a library of 100,000 
model SEDs in 7 bands (FUV, NUV, u, g, r, i, z,) by generating 
model galaxies with a range of ages, star formation histories, metal- 
licities and dust attenuation strengths. For a given galaxy, they eval- 
uated the;^'^ goodness-of-fit by comparing the observed photometry 
with each model, determined the relative weight for each model, 
and finally built a probability distribution function (PDF) for the 
model parameters. The average of the PDF is adopted as the a nom- 
inal estimate of the parameters. 

S07 show that the derived dust attenuation (the optical depth 
Tv) is sensitive to the assumed prior distribution of Tv in the model 
library, so the prior distribution should be as realistic as possible. 
To solve this problem, we estimate approximate values of the at- 
tenuation Tv(gas) fro m the B aimer decrem ent F(Ha')/F(H/J) within 
the SDSS fiber and the lCalzetti eTai] j2000l) extinction curve, when 
both F(Hq') and F(H/3) have S/N>3. When either F(Hq-) or F(H/3) 
has 5/^ < 3, we adopt Tyistar, fiber)/0.44 from the measure- 
ments of the attenuation of the stellar continuum of the galaxy 
provided in the MPA/JHU catalog (tauv^cont). This parameter is 
obtained by fitting Bruzual & Chariot (2003) population synthesis 
models to the stellar continuum; the reddening may then be esti- 
mated by determining the extra "tilt" that must be applied to the 
models in order to fit the shape of the observed spectrum. 

The prior distribution of Tv(gas) and fj (Tv(star) = Tv(gas) x 
/j) values we adopt in our model library is tuned to reproduce 
the Tyistar, fiber) predictions for the whole HI sample. We adopt 
for Ty(gas) a Gaussian distribution peaked at 1.78, with a width 
cr = 0.55, and we ad opt for /i a Gaussian ditribution peaked at 0.44 
jCalzetti et alj2000l) . with a width cr = 0.4. The Gaussian distribu- 
tion for Ty(gas) is trimmed so that it only spans values between 
and 4, and the Gaussian distribution for // is trimmed so that it only 
spans values between and 1. Here we assume that the Ty(star) 
does not change significantly throughout the galaxy. 

In reality, the dust distribution in galaxies is more complicated 
than our simple model assumes. One issue is that for a typical 
galaxy, the relative attenuation between t he inner disk, out er disk 
and bulge depends strongly on inclination jTufi's et alj2004l) . How- 
ever, this effect becomes most notable when the inclination is large 
(axis ratio b/a < 0.4), so this problem is much mitigated by ex- 
cluding galaxies with b/a < 0.4 from the samples (Section l2.lt 
( I Yip et al.ll2O10l) . In the absence of more detailed information on 
the dust distribution in each galaxy, we believe our method gives 
results that are as accurate as possible. 

To test the robustness of the sSFR derived from our SED fit- 
ting, we also measured sSFR directly from the SDSS fiber spec- 
trum, and we compared this to the sSFR estimated within the 
3 " aperture using our SED fitting technique. The SFR derived 
from the spectrum is calculated from the Hey luminosity, and cor- 
rected for dust using the Balmer decrement F(Ha)/F(Hj8) and the 



ICalzetti et all 1 I2OOOI) extinction curve. The stellar mass inside the 
SDSS fiber is taken from the MPA/JHU catalog. For SED fit- 
ting, we take the SDSS fluxes within the fiber from the SDSS 
archive (fibercounts), measure the GALEX NUV fluxes and con- 
volved SDSS u band and GALEX FUV fluxes within 6 " around 
the galaxy (approximately the FWHM of the GALEX NUV PSF). 
Then we use the u band fluxes (within the fiber aperture and from 
the convolved images) to normalize the GALEX FUV and NUV 
fluxes, so that we get consistent SED from the 7 bands within the 
3" aperture. 

We see from Figure [7] that the two estimates generally agree 
with each other, with a Icr scatter of ~0.39 dex. The offset between 
the two estimates exhibits weak systematic trends with both log 
SFR and Ty. The most discrepant results are obtained for galaxies 
with very low specific star formation rates, where the UV luminosi- 
ties are low and may trace populations of stars that are not properly 
accounted for in standard population synthesis models (Com'oy & 
Gunn 2010). We have also tested the effect of varying the adopted 
priors for Ty . We confirm that the prior that is most similar to the 
real distribution (as calculated from the spectrum), gives the best 
fit. 



3.4 Morphological parameters 

We closely follow the procedure described in iLotz et al. I ( l2004h . 
hereafter L04, to derive the asymmetry (A) and smoothness (5 ) pa- 
rameters from the SDSS g-band images. Both A and S are mea- 
sured inside an aperture equal to 1.5 times the petrosian radius (r,,). 
All neighboring objects are masked. 

A is a measure of the difference between a given galaxy image 
and the image rotated by 180 degrees about the object center (the 
center is determined by minimizing A). A higher value of A means 
that the galaxy is more asymmetric. 5 is a measure of the difference 
between a given galaxy image and the image smoothed with a 0.2r,, 
wide boxcar kernel. A higher value of 5 implies that the galaxy has 
a more clumpy morphology on scales equal to the kernel size. A and 
S are also calculated using regions of blank sky in the vicinity of the 
galaxy; the average values of A and S calculated for the background 
are subtracted from the measurements we make for the galaxy. 

To make sure that our measurements are consistent with L04, 
we measure A and 5 for 41 galaxies from L04 with images avail- 
able from the SDSS. Our results are shown in Figure [8] We see 
that our parameters and those from L04 follow similar trends along 
the Hubble sequence. There is quite close agreement between the 
A values, but our S values are typically four times smaller than 
those of Lotz - this is likely to be a reflection of the different qual- 
ity/resolution of the images or different details in the image pro- 
cessing steps in the two cases. Nevertheless, our test shows that the 
relative trends in 5 are the same for both sets of measurements, 
so these parameters are a useful diagnostic of relative changes in 
galaxy morphology. 

We visually check all the images in our sample using the 
SDSS Imaging Finding Chart TooljJ and we identify galaxies that 
are clearly undergoing a merger event. We find no enhancement in 
the merger fraction for the galaxies in the HI sample (9, 8, and 7 
galaxies out of a total of 519 in the HI, Cm,, and CM,j,uvr samples 
exhibit clear signs of a merger or interaction in our images). 



' http://cas.sdss.org/astro/en/tools/chart/chart.asp 
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4 RESULTS 



As shown in Figures 2 and 3, galaxies from the HI sample are bluer, 
and have larger sizes and later type morphologies compared to C^. 
galaxies. This is not surprising, because it is well known that galax- 
ies that are more actively star-forming also contain more gas. The 
CM,,Nuv-r and CM,.Nuv-r.fit galaxies, on the other hand, are matched 
to the HI sample both in stellar mass and in global NUV-r colour, 
so if HI gas fraction and star formation activity track each other 
very closely, one might expect the galaxies in these two samples to 
have identical properties. In the following sections, we will show 
that this is not true. We divide the HI sample into four stellar mass 
bins, and study how sizes and colour gradients vary with atomic 
gas mass fraction. We also compare properties of the correspond- 
ing control galaxies along the same sequence. The main purpose of 
comparing our results with those derived from control samples is to 
isolate those trends that can be attributed to increasing HI content, 
rather than to any other correlated property, such as stellar mass 
or global star formation rate. Because the samples have all been 
matched in redshift, the noise in the measurements of luminosity 
and colour necessary to estimate quantities such as stellar mass, 
will be identical in all the comparison samples. 



4.1 HI mass fraction and the optical sizes of galaxies 

Figures[3]and|4]demonstrate that the galaxies in the HI sample have 
larger average sizes than both Cm, and Cm-tMuv-r control sample 
galaxies. The mean value of Rio(i) of the HI sample is larger by 1.2 
and 0.3 kpc than the mean values of Rso(i) of the Cm, and CM:Nuv-r 
control sample galaxies respectively. In Figure |9l we plot the rela- 
tions between the half-light radius i?5o and M(H1)/M, in four dif- 
ferent bins of stellar mass for galaxies in the HI sample. The mean 
values of M(H1)/M, for each stellar mass bin are marked as crosses 
of different sizes at the bottom of the plot. These have been derived 
by Catinella et al (2010) using GASS survey galaxies for which HI 
masses have been measured down to a limiting M(H1)/M, limit of 
~ 1.5%. The majority of our Hl-detected galaxies have atomic gas 
fractions that are above the average value - this is not surprising 
since galaxies detected in the ALFALFA survey make up the bulk 
of our sample. 

We see that the half-light radius R50 measured in the g-band 
increases as a function of M(H1)/M, at a given value of M,. This 
result is consistent with the scaling relations published in Zhang 
et al (2009) and CIO, which clearly showed that M(H1)/M. corre- 
lates strongly with stellar surface density /j, (i.e. with galaxy size 
at a fixed value of the stellar mass). Indeed, both studies find that 
M(H1)/M, can be best predicted using a combination of colour and 
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jj,. As discussed in Zhang et al (2009), these scahngs can be un- 
derstood in term s of the Kennicutt-Schmidt law of star formation 
jKennicuttll 199 Sh . which states that the surface density of star for- 
mation scales with the surface density of cold gas as a power-law 
with slope ~ 1.4. A star formation law of this form leads to the 
expectation that 

log M(H I) /M* = alogfi, + blogS FR/M, + c, (1) 

where a, b and c are constants. The NUV-r colour is an excellent 
proxy for SFR/M, (especially in the blue sequence), so this leads 
to a prediction of a linear relation linking HI mass fraction, NUV-r 
colour and stellar surface mass density. The two right hand panels 
of Figure |9] show that sizes of CM,,Nuv-r and CM,,Nuv-r4u control 
galaxies exhibit the same increase as for the Hl-detected galaxies. 
This supports our hypothesis that the scaling of galaxy size with HI 
mass fraction arises as a consequence of the star formation rate-gas 
surface density relation. 

In Figure [Tol we analyze trends in the ratio of the ^-band and 
i-band half-light radii as a function of HI mass fraction. R50 de- 
fined in the ^-band is always slightly larger than R50 defined in 
the /-band, because the light from younger stellar populations is 
generally spread over a larger effective radius than the light from 
older stellar populations. The ratio between the g and /-band radii 
increases as a function of gas fraction, and is also larger for more 
massive galaxies at a fixed value of M(HI)/M.. As we will show in 
the next section, this difference in g- and /-band half-light radii is 
also found when we analyze g - / colour gradients. In the panels 
to the right, we plot size differences for the three control samples. 
Neither the Cm, nor the CM,.Nuv-r control samples show similar 
trends to the HI sample. The size differences do increase for more 
massive galaxies in the Cm,.nuv-i;ii, control sample, but the effects 
are weaker than those observed for the sample with real HI mass 
measurements. 

One might worry whether these trends arise because Hl-rich 
galaxies are more dusty. We analyzed the correlation between op- 
tical depth Ty and HI mass fraction, and find that Ty does not in- 
crease with M(HI)/M,. In addition, control galaxies matched in 
stellar mass and colour do not have smaller Ty values than galaxies 
from the HI sample. We are thus led to the conclusion that there is 
an intrinsic correlation between M(HI)/M, and R50(g)/R50(i). 



4.2 HI mass fraction and (g - /) colour gradients 

In the top two panels of Figure (TT] we analyze trends in the inner 
and outer g - i colours of the galaxies in our four samples. In the 
bottom panel, we look at trends in (g - /) colour gradients. Our 
results are as follows: 

Inner colours: The inner g - / colour of a galaxy depends 
strongly on stellar mass, in the sense that more massive galaxies 
have redder inner colours. The inner g - / colour depends only 
weakly on HI mass fraction at a given value of M,. We see very 
similar trends in inner g - / colour for control samples matched in 
NUV - r colours. We conclude that the inner colour of a galaxy is 
primarily sensitive to its stellar mass. 

Outer colours: The behaviour of outer g - / colour as a func- 
tion of stellar mass and HI fraction is very different to that of the 
inner colour. It depends strongly on HI mass fraction, in the sense 
that more Hl-rich galaxies have bluer outer colours. At a fixed value 
of the HI mass fraction, there is no dependence on stellar mass. We 
see very similar trends in the CM,,Nuv-r and CM.j^uv-r./j, control 
samples, but not in the Cm. control sample. This indicates that the 



outer colour of a galaxy is primarily sensitive to how much ongoing 
star formation there is in the galaxy. 

Colour gradients: The g - / colour gradient of a galaxy de- 
pends strongly on its HI mass fraction, in the sense that more Hl- 
rich galaxies become bluer on the outside relative to the inside. At 
a fixed value of M(HI)/M., more massive galaxies have stronger 
colour gradients. None of the control samples exhibits exactly the 
same trends as the HI sample. In the CM,.Nuv-r and CM„Nuv-r,ii. 
samples, the colour gradients for massive galaxies behave in much 
the same way as for the sample with HI measurements. For low 
mass galaxies, however, there is no trend in colour gradient with 
log M(HI)/M, for the control galaxies. Our main conclusion, there- 
fore, is that there is an intrinsic link between higher HI fractions 
and bluer outer disks. 

4.3 HI mass fraction and gradients 'in{NUV - r) colour and 
sSFR 

In this section we analyze trends in the inner and outer NUV-r 
colours and colour gradients as a function of HI mass fraction. The 
analysis is similar to that of the g - i colours in the previous sec- 
tion, except that in this case, the average inner and outer colours are 
computed from stacked images due to the lower S IN of the NUV 
images. As we will show, the main advantage of studying NUV-r 
colours is the fact that the UV is a much more sensitive probe of the 
presence of young stars than the g-band, so that both the trends as a 
function of M» and M(HI)/M,, and the differences between the HI 
sample and the control samples found in the previous section are 
much enhanced. 

The results shown in Figure[T2]confirm essentially all the con- 
clusions discussed in the previous section. The main qualitative 
difference we find by studying NUV-r colours rather than g - / 
colours is that the behaviour of the colour gradients and the outer 
NUV-r colours as a function of HI mass fraction now differ more 
significantly from those exhibited by both the CM*.Nuv-r and the 
CM,.Nuv-r,^i, control galaxies. This strengthens our conclusion that 
the trend in the colour of the outer disk is fundamentally linked to 
the increase in the HI content of the galaxy. 

We note that Figures II Hand [T2]both show that it is the most 
massive galaxies with stellar masses ~ 10" Mq, and with relatively 
low HI gas fractions of ~ 10%, that have the strongest colour gradi- 
ents. Low mass galaxies have the highest HI gas fractions, but their 
colour gradients are relatively weak. It is interesting to investigate 
the morphology of the UV light for both kinds of systems. Fig- 
ure[T3]shows optical and NUV images for the 18 Hl-sample galax- 
ies with the most negative A„_j{NUV - r) < -1.1 as well as the 18 
Hl-sample galaxies which have the highest values of M(HI)/M, in 
the HI sample (we point out that, because the HI sample is stellar 
mass limited at 10'" Mq, we are missing a population with the high- 
est HI mass fractions {M(HI)/M,) at a fixed M(HI) ). The galaxies 
with the most negative A^_j(NUV - r) all have very red central re- 
gions and a blue outer disk. Many of them have bars, and the blue 
light tends to be located in a UV-bright ring around the bar (up- 
per panels of Figure|13l>. Galaxies with the highest HI gas fractions 
are usually blue throughout the whole galaxy (lower panels of Fig- 
urell3b. 

Finally, Figure [14] shows plots equivalent to Figures 1111 and 
ll2lfor the specific star formation rates we derive from fitting models 
to the UV through optical SEDs derived from the stacked images 
(see section 3.3). We caution that we have to make assumptions 
about the likely range of dust attenuation in our galaxy sample, 
as well as assumptions about how the dust might be distributed 
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Figure 9. Relation between galaxy size (from the g band) and HI mass fraction for galaxies in the HI sample (blue circles), Cm* control sample (red diamonds), 
C-Mt.NUV-r control sample (green triangles) and CMt,nuvr,^t control sample (cyan squares). The control galaxies have no HI mass measurements, and are just 
plotted in the same M(HI)/M, bin as the corresponding galaxies in the HI sample. Symbols of different sizes are for different stellar masses as indicated at the 
corners of the plots. The error bars show the r.m.s. deviation in the mean galaxy size calculated through bootstrapping. The crosses at the bottom of the plot 
mark the location of the mean gas fraction for galaxies in each stellar mass bin as given in CIO. 
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Figure 10. As in Figure 8, except for the ratio between the g and ('-band half-light radii. 
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Figure 11. Relation between inner g - i colour, outer g - i colour and A„_,(g - i) and HI mass fraction for the HI sample (left column, plotted as blue circles). 
Cm* control galaxies (second column, plotted as red diamonds) , the CM».jv(/v-r control galaxies (third column, plotted as green triangles) and the Cut.NUV-r.^* 
control galaxies (fourth column, plotted as cyan squares). The control galaxies have no HI mass measurements, and are plotted in the same M(HI)/M» bin as 
the corresponding galaxies in the HI sample. The sample is divided into 4 stellar mass bins, and plotted with symbols of different sizes as indicated in Figure|9] 
The error bars show the r.m.s. deviation in the mean values calculated through bootstrapping. The crosses denote the mean gas fraction for galaxies in each 
stellar mass bin (from CIO). 
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Figure 12. Relation between inner NUV — r, outer NUV — r, ha-i(NUV — r) and HI mass fraction for the galaxies from the HI (blue circles) sample, and Cm, 
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throughout the galaxy. We thus regard these results as indicative, 
rather than definitive. The main point we wish to make with this 
plot is that with plausible assumptions about extinction, we find 
that outer disks are younger, with a higher ratio of present-to-past 
averaged star formation in galaxies with higher HI mass fractions. 
The same effect is not seen as clearly or strongly in the control 
samples. This implies that the outer disks of galaxies with higher HI 
mass fractions are are growing more rapidly, i.e. that these galaxies 
are forming from the "inside-out". 

We have estimated the uncertainty in our estimated sSFR aris- 
ing from the age-metallicity degeneracy. The offset between the 
estimated value of the sSFR will be less than 0.2 dex for galaxies 
with half-solar and twice solar metallicities. Note also, that all the 
galaxies in this analysis hav e M, > IO'^Mq; the mass-metallicity 
relation taallazzi et al.l2005h implies that there will be no very low 
metallicity galaxies included in our analysis. 



We conclude that the correlation of galaxy asymmetry and 
smoothness with HI mass fraction arises because galaxies with 
higher HI gas fractions have younger stellar populations. Reichard 
et al (2009) show that there is a strong link between galaxy asym- 
metry (in their case, as measured by the m = 1 azimuthal Fourier 
mode) and the age of its stellar population. They found that this link 
was independent of the other structural properties of the galaxy. 
Our analysis of the control samples appears to indicate that this 
link is also independent of HI content. Lopsidedness/asymmetry 
is believed to trace a non-equilibrium dynamical state caused by 
mergers, tidal interactions, asymmetric accretion of gas, or a sym- 
metries related to the dark matter halo teoumard etXll2005l) . As 
discussed in section 3.4, we found no enhancement in the fraction 
of galaxies undergoing mergers in the HI sample relative to control 
samples, so it is perhaps not too surprising that we find no intrinsic 
correlation between lopsidedness/asymmetry and HI content. 



4.4 HI gas fraction and higher order morphological 
parameters 

FigurefTSlshows A and 5 as a function of HI fraction for the galax- 
ies in the HI sample (blue), and for galaxies in the Cm, (red), 
CM,,Nuv-r (green) and CM,.Nuv-r.fi, (cyan) control samples. Galax- 
ies in the HI sample are significantly more asymmetric and less 
smooth than galaxies in the Cm, sample. There is no significant 
difference between the HI, CM,,Nuv-r and CM,,Nuv-r,ii, samples. In 
Figure [T6l we plot A and S as a function of HI mass fraction for 
galaxies in the HI sample (blue) and in the CMt.Nuv-r control sam- 
ple (green) for different stellar mass bins. We see that at fixed stellar 
mass, A and 5 both increase as a function of M(HI)/M» for galaxies 
in the HI sample. At fixed stellar mass, galaxies in Cm^nuv-v sam- 
ple have similar values of the asymetry parameter A and smooth- 
ness parameter S (except at the very highest HI mass fractions). 



5 DISCUSSION AND CONCLUSIONS 

We have analyzed a sample of galaxies from the parent sample 
of the GASS survey, for which HI mass measurements are avail- 
able from the GASS survey itself and from ALFALFA. At a given 
value of M,, our sample consists primarily of galaxies with larger- 
than-average values of M(HI)/M, . We use a combination of optical 
photometry from the SDSS and FUV,NUV images from GALEX to 
study properties of the 2-zone stellar populations of these Hl-rich 
systems, and we investigate how they vary with HI mass fraction. 

We constructed three control samples for comparison with the 
galaxies with HI mass measurements: one sample is matched in 
stellar mass and redshift (Cm*), the second one is matched in stellar 
mass, NUV - r colour and redshift {CM:Nuv-r) and the third one 
is matched in stellar mass, NW - r colour, surface mass density 
and redshift (Cm^j^vv-A- We then generated self-consistent 7-band 
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Figure 13. SDSS and GALEX NUV images of galaxies from tlie HI sample which have the most negative (NUV - r) colour differences (Ao-i(NUV — r)\-l.l) 
(the first 3 rows), and of galaxies that have the highest HI mass fractions (M(HI)/M,) (~ 1) (last 3 rows). For each image pair, the left image is the SDSS g, r, i 
colour composite image, and the right image is the NUV map for the same galaxy. 



photometry (FUV, NUV, ii,g,r,i and z) for the galaxies in all three 
samples and derived inner colours, outer colours, asymmetry and 
smoothness parameters for all the galaxies. We also generated four 
sets of stacked images in different bins of stellar mass and HI gas 
fraction, and used these to derive inner and outer specific star for- 
mation rates using an SED-fitting technique. 

The main purpose of comparing our results with those derived 
from control samples is to isolate those trends that can be attributed 
to increasing HI content, rather than to any other correlated prop- 
erty, such as stellar mass or global star formation rate. As expected, 
Hl-rich galaxies differ in nearly all their properties from galaxies of 
the same stellar mass selected without regard to HI content. Most 
of these differences are attributable to the fact that galaxies with 
more gas are bluer and more actively star-forming. In order to iden- 
tify those galaxy properties that are truly causally connected with 
HI content, we also compare our results derived for the HI sample 
with those derived for galaxies with the same distribution of both 
stellar mass and NUV-r colour. Finally, Catinella et al (2010) have 
demonstrated that the best predictor of the HI gas fraction is a lin- 



ear combination of NUV-r colour and the logarithm of the stellar 
surface mass density of the galaxy. At a fixed colour, galaxies with 
lower stellar mass densities have higher HI gas fractions. This mo- 
tivates comparison with a control sample that is matched in stellar 
mass, NUV-r colour and stellar surface mass density. 

We find that the only photometric property that is clearly at- 
tributable to increasing HI content and not to the total amount of 
star formation in the galaxy, is the colour gradient of the galaxy. 
There appears to be a causal link between HI fraction and outer 
disks that are bluer with respect to the inner part of the galaxy. The 
same effect also manifests itself as a systematic increase in the ratio 
of size of the galaxy measured from SDSS images in the g-band to 
that measured in the i-band. The colour gradients are significantly 
stronger if one analyzes UV/optical colours than if one uses g - i 
as a tracer. Our SED modelling demonstrates that as the HI frac- 
tions of galaxies increase, the specific star formation rates in their 
outer regions become higher with respect to their inner regions. 
This means that the outer regions of Hl-rich galaxies are younger. 

Control samples that are matched in size in addition to colour 
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Figure 14. Relation between inner sS FR, outer sSFR, Ao_,(.v5fS) and HI mass fraction for galaxies from the HI (blue circles) sample, and from the Cu, 
(red diamonds), CM».mnT (green triangles) and CM»./imT,fi« (cyan squares) control samples. The sSFR are derived from SED fitting as explained in the text, and 
the SEDs are measured from the stacked images. Note that the control galaxies have no HI mass measurements, and are plotted in the same M(HI)/M, bin as 
the corresponding galaxies in the HI sample. Symbols of different sizes indicate different stellar masses as denoted in Figure|9](larger symbol corresponds to 
higher stellar masses). The crosses denote the mean gas fraction for galaxies in each stellar mass bin (from CIO). 
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Figure 15. Relation between the median values of the morphological parameters asymmetry (A) and smoothness (5 ) and HI mass fraction for galaxies in 
the HI sample (blue symbols and lines). The median A and 5 values of the corresponding Cm* galaxies are plotted in red; results for CM*,NUV-r galaxies are 
shown in green; results for CMt,NVV-r,ii* galaxies are shown in cyan. A median value is plotted only when there are more than 6 galaxies in that bin. Error bars 
show the r.m.s. deviation calculated through bootstrapping. 



and stellar mass, tend to exhibit stronger colour gradients, but the 
observed scalings are still weaker than for the sample with real HI 
mass measurements. This indicates that colour gradients and size 
are partially co-variant, as might be expected if infalling gas was 
recently added to the outside of the disk in the most Hl-rich galaxies 
in our sample. 

The "inside-out" picture of disk galaxy formation has com- 
monly served as a basis for semi-analytic modelling the chemical 
evolution of disk galaxies (e.g. Chiappini et al. 1997), as well as 
the formation of disk galaxies in the context of Cold Dark Matter 
cosmologies (e.g. Kauffmann 1996; Van den Bosch 1998; Boissier 
& Prantzos 1999; Dutton 2009). In this scheme, gas in the halo sur- 
rounding the galaxy cools, falls onto the galaxy while conserving 
its angular momentum, and fuels star formation in the disk. Gas ac- 



creted at late time has higher specific angular momentum, and so 
settles in the outer region of the galaxy. However, hydrodynamical 
simulations of the formation of disk galaxies indicate that angular 
momentum is not conserved and disks do not always form from the 
inside out (Sommer-Larsen et al 2003; Robertson et al 2004). 

Direct observational evidence that disk galaxies form from the 
inside out has thus far been quite difficult to establish. Pohlen & 
Trujillo (2006) found that 60% of local disk galaxies have inner ex- 
ponential profiles followed by a steeper (downward bending) outer 
exponential profile, while 30% have a shallower (upward bending) 
outer profiles. Azzollini et al (2008) find that for a given stellar 
mass, the radial position of this break has increased with cosmic 
time by a factor of 1.3-I-/-0. 1 between z = I and 0, in agreement 
with a moderate inside-out growth of disk galaxies over the last ~ 8 
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Figure 16. Relation between A and 5 and HI mass fraction for galaxies from the HI (blue circles) sample and from the Cut^NUV-r (green tiiangles) control 
sample. The sample is divided into 4 stellar mass bins, indicated by symbols of different sizes as described in the corner of each panel. The error bars show the 
r.m.s. deviations in the mean values of A and S calculated through bootstrapping. The crosses denote the mean gas fraction for galaxies in each stellar mass 
bin (from CIO). 



Gyr. Cresci et al (2009) find evidence for evolution in the zero point 
of the TuUy-Fisher relation of galaxy disks at z ~ 2 with respect to 
local samples that is in agreement with cosmological simulations 
of disk galaxy formation. 

Coherent studies of the metallicity, colour and star formation 
rate gradients in local galaxies have been few and far between. 
Zaritsky, Kennicutt & Huchra (1994) performed a careful study 
of abundance gradients derived from HII regions in a sample of 
53 nearby spiral galaxies. They found that the slopes of the ra- 
dial abundance gradients when rescaled to units of dex/isophotal 
radius, do not correlate with the luminosity or Hubble type of the 
galaxy. They did find evidence that galaxies with bars have flat- 
ter abundance gradients, which suggests that gas inflows may be 
import ant in understan ding the time evolution of metallicity gra- 
dients. Ide JongI ( 1 19961) investigated the effects of star formation 
history and the metallicity on the colour profiles of 86 face-on spi- 
ral galaxies, and found that the outer regions of disks are on aver- 
age younger and of lower metallicity. Munoz-Mateos et al (2007) 
studied specific star formation rate profiles of a sample of 160 
nearby spiral galaxies from the GALEX Atlas of Nearby galaxies 
jGil de Paz et al.ll2007l) . They found a large dispersion in the slope 
of the sSFR profiles, with a slightly positive mean value, which 
they interpreted as implying moderate net inside-out disk forma- 
tion. They found that in a minority of galaxies, the scale length 
of the UV light is significantly larger than the scale length of the 
near-lR light, and interpreted these galaxies as having undergone 
an episode of enhanced recent growth. 

In this paper we have uncovered a direct link between the HI 
content of a galaxy and its specific star formation rate profile. This 
lends further credence to the idea that galaxies where the outer disk 
is very blue and actively star-forming compared to the inner disk 
have recently accreted gas. The fact that there is no intrinsic cor- 
relation between the HI fraction and the measured asymmetry of 
the optical light of the galaxy, suggests that the gas was accreted 
smoothly and not in discrete units. What is still lacking, however, 
is clear proof that the gas in Hl-rich galaxies has been added to 
the outer disk recently. It will also be extremely important to ob- 
tain data that will allow us to link the cold gas content, specific star 
formation rates and metallicity gradients of disk galaxies in a self- 
consistent way. Only then, may we hope to fully constrain models 
of disk galaxy formation and chemical evolution. 



We are currently engaged in an observational program to ob- 
tain accurate and homogeneous molecular gas masses for a subset 
of 300 galaxies from the GASS sample (Saintonge et al, in prepara- 
tion; see |http://www.mpa-garching.mpg.de/COLD~GASS/| l. In ad- 
dition, we are obtaining long-slit spectra for the same galaxies us- 
ing the MMT telescope (the 6.5m MMT telescope on Mt. Hop- 
kins, AZ). Early results indicate that specific star formation rate 
and me tallicity gradients in Hl-rich galaxies are indeed tightly cor- 
related l lMoran et alteoid) . We intend to quantify this in more de- 
tail with larger samples of galaxies in our future work. 
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APPENDIX A: A TEST OF THE ROBUSTNESS OF OUR 
2-ZONE UV/OPTICAL COLOUR MEASUREMENTS 

In this paper, UV/optical colours are measured by convolving the 
SDSS images to match the poorer resolution of the GALEX NUV 
images. In some of our galaxies, the half-light semi-minor axis Z750 
becomes similar to or smaller than the size of the ~ 4 pixel PSR 
of the GALEX NUV images. In this regime, the 2-zone colour 
and SFR measurements may no longer be valid. In addition, the 
GALEX MIS images are relatively shallow, so the 2-zone mea- 
surements of the NUV fluxes may have large errors if the signal-to- 
noise is low. 

In order to quantify these effects more precisely, we have per- 
formed simulations by transforming the images of nearby galaxies 
so that they have similar magnitudes and apparent size ranges to the 
galaxies in the samples studied in this paper. We s elect 52 galaxies 
from t he GALEX Nearby Galaxies Survey (NGS. lGil de Paz et al.l 
j2007) ) with images available from the SDSS. We only select iso- 
lated galaxies with apparent sizes (D25) less than 5 arcmin, where 
the image is not contaminated by foreground or background ob- 
jects. 

We begin by measuring NUV - r and Ao-i(NUV - r) from 
these images (shown in purple in Figure A3). The nearby sample 
turns out to have a similar distribution of these two parameters as 
the GASS parent sample (shown in black on in Figure A3). 

We then create a library of 60,000 simulated images. For each 
nearby galaxy image, we select a value of b<,o and r-band appar- 
ent magnitude at random from the distributions shown in black in 
Figure A3. We also select a NUV image exposure time at random, 
again based on the real distribution of exposure times for the GASS 
parent sample. The NUV and r band fluxes are always rescaled 
so that the total NUV - r colour is conserved. After rebinning in 
size, the NUV and r band images are convolved with the GALEX 
and SDSS PSFs and background and poisson noise are added. Fig- 
ures |AT] and |A2] show two examples of the resulting simulated im- 
ages. Figure lATI shows a typical galaxy, which is blue on the outside 
(with negative A„_,(A't/K-r)), while Figure lA2l shows a galaxy that 
is blue on the inside (with positive A„_,(A't/V - r)). 

After creating our library of 60,000 simulated images, we 
measure photometric parameters following the same steps de- 
scribed in Section im The distribution of output NUV-r colours, 
r-band apparent magnitudes, semi-minor axis size (fojo) and NUV- 
r colour gradient values are plotted as red histograms in Figure A3. 
We can see that the output values do differ from the input values. 

To quantify this in more detail, the left panel of Figure IA4I 
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shows how the difference between output and input colour gradi- 
ents change as a function of for galaxies with high S IN UV 
images {cr{NUV - r) < 0.1). When bso > 4 pixels (which corre- 
sponds to 6 "), the output Ac_j(NUV - r) is quite close to the input 
value (systematically high by ~0.05 mag on average). The scatter 
in the recovered gradient is less than 0.2 mag. When ^150 < 4 pixels, 
the mean value of diff(A„_j(NUV-r)) rises steeply and the scatter 
in the recovered value also increases. We thus impose a lower limit 
on of 4 pixels when we analyze UV/optical photometry in this 
paper (see Sect. 3.2). 

The right panel of figure lA4l shows the difference between the 
output and input colour gradients as a function of the errors on 
the NUV - r colour, for galaxies with fcjo > 4 pixels. The average 
offset remains constant at ~0.05 mag when -log criNW -r) > 0.5, 
and then rises steeply towards lower values of -log a{NUV — r). 
However, when -log cr{NUV - r) < I (or cr(NUV - r) > 0.1), the 
scatter in the recovered gradient is already about 0.2 mag. This is 
because the outer UV light has lower surface brightness than the 
inner light. We conclude that o-(NUV - r) < 0.1 mag should be 
adopted as a lower 5 /N limit for robustly measuring total and 2- 
zone NUV-r colours. 
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Figure Al. Original and simulated images of the galaxy NGC 5701. The left-top panel shows the original image, the top-right panel shows a simulated image 
with high S IN, the bottom-left panel shows a simulated image with low S IN, and the bottom-right panel shows an unresolved simulated image with high 5 IN. 
The two ellipses on each galaxy image mark the inner and outer regions of the galaxy. 




Figure A2. As in the previous figure, except for galaxy NGC 4314, which has a blue core. 
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Figure A3. Distributions of NUV - r colours, r band appai'ent magnitudes, semi-minor axis radii (fojo) measured in the r-band, and colour differences 
(Ao_,(A'f/y — r)) for the GASS parent sample (black), for the original NGS galaxies (blue), and for the simulated images (red). 
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Figure A4. The average difference between the A„_,(A'J/V - r) colour for the simulated galaxy images and the original galaxy images is plotted as a function 
of half-light radius of the galaxy in arcsec (left), for cases where it(NUV — r) < 0.1. In the right panel, the average difference is plotted as a function of 
NUV - r colour errors, for cases where fejo > 4 pixels. The solid lines show the average difference in each b^Q bin, the dashed line shows the median, and the 
error bars show the r.m.s. scatter in the difference around its mean value. The dotted lines mark where the offset equals 0, and the dot-dashed lines mark the 
location of our recommended size and error thresholds discussed in Sect. 3.2. 



